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Neutron diffraction, magnetization, resistivity, and heat capacity measurements on the 6H- 
perovskite Ba3BiRu209 reveal simultaneous magnetic and structural dimerization driven by strong 
magnetoelastic coupling. An isostructural but strongly displacive first-order transition on cooling 
through T* — 176 K is associated with a change in the nature of direct Ru— Ru bonds within RU2O9 
face-sharing octahedra. Above T* , Ba3BiRu20g is an S = 1 magnetic dimer system with intradimer 
exchange interactions Jo/fce ~ 320 K and interdimer exchange interactions J' /ks = —160 K. Below 
r*, a spin-gapped state emerges with A ~ 220 K. Ab initio calculations confirm antiferromagnetic 
exchange within dimers, but the transition is not accompanied by long range-magnetic order. 

PACS numbers: 61.05.fm, 75.10.Pq, 75.30.Et, 75.30.Kz, 75.40.-s, 75.47.Lx 
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Quantum cooperative phenomena involving charge, 
spin, orbital, and lattice order parameters are among the 
most explored areas of modern solid-state physics. A 
particularly rich vein of quantum cooperative phenom- 
ena are low-dimensional magnetic systems featuring mo- 
tifs such as dimers, chains, ladders, and plaquettes, where 
coupling between spin and lattice degrees of freedom can 
give rise to magnetoelastic effects such as the spin-Peierls 
transition in S" = 1/2 Ising-chain antiferromagnets [ij. 
The most thoroughly studied are based on = 1/2 3d 
transition metal cations such as Cu^+, and Ti'^+. 
The less common S = 1 systems remain relatively ne- 
glected, although S — 1 ruthenates (low-spin Ru''+), in 
particular, show fascinating properties including super- 
conductivity (Sr2Ru04[3|), non- Fermi liquid behaviour 
(La4Ru60i9^3]) and low-dimensional character (the fa- 
mous spin-gapped Haldane phase Tl2Ru207[3|)- 

This report concerns Ba3BiRu209, first reported by 
Darriet et al.[§\ as a 6H-perovskite with a small mon- 
oclinic (C2/c) distortion. It contains BiOg octahe- 
dra sharing vertices with Ru2 0g face-sharing octahedral 
dimers, and Ba^+ ions occupying high-coordinate A sites 
(Fig. [1]). It forms part of the Ba3i?Ru20| series, where 
i? is a rare-earth or 3d transition metal G,?], indium's'] 
or zirconiumj^. The rare-earth R cations usually have a 
3-1- oxidation state, giving Ruj"*" dimers; but Tb, Pr, and 
Ce have 4-1- oxidation states, giving Ruj"*" dimers, and 
shrin king the unit cell due to the reduced ionic radius 
of R [10| (Fig. [IJ . Our experimental lattice parameters 
for Ba3BiRu209 strongly suggest that it belongs in this 
Ruj^ category. 

Polycrystalline Ba3BiRu20g was synthesized by solid- 
state reaction. Neutron powder diffraction (NPD) data 
were collected on the instrument Wombat at OPAL, 
ANSTO, from 2 - 300 K using 2.9609 A neutrons; and 
at 2 K on the HRPD diffractometer at ISIS. Structure 
refinements via the Rietveld method were carried out us- 
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FIG. 1. (Color online) Ionic radu ^ of J?^+ {R = Bi, Y, 
rare-earth) in 6-fold coordination vs. unit cell lengths a, b/y/S 
(filled circles) and c (open squares) for Ba3i?Ru209 at 298 K. 
Error bars are smaller than symbols. Adapted from Doi et al. 

and references therein, with i? = Bi from the present work. 
The refined (HRPD@ISIS) 2 K structure of BasBiRuaOg is 
shown with BiOe octahedra (violet), RU2O9 dimers (grey), 
Ba?^ (green, and O^^ ions (red). Polyhedral units are ex- 
panded and drawn with 90% probability thermal ellipsoids. 



ing GSASfTT'l with the EXPGUI front-end Pi . Magnetic 
susceptibility, electrical resistivity and heat capacity were 
measured in a Quantum Design PPMS. Ab initio calcu- 
lations were performed in the generalized gradient ap- 
proximation (GGA) using the Vienna ab initio Simula- 
tions Package (VASP 5.2) [13]. A supercell containing two 
primitive cells (60 atoms) was used with standard PAW 
potentials [3], a fc-mesh of 162 points in the irreducible 
Brillouin zone wedge and a cutoff energy of 450 eV. Total 
energy converged to within 10~^ eV. 

Fig. [2] displays the temperature dependencies of the 
unit cell volume, lattice parameters, and intradimer 
^Ru-Ru distance. The volume shows normal thermal con- 
traction down to ~180 K, then increases rapidly, mostly 
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FIG. 2. (Color online) Temperature dependencies of a) unit 
cell volume V, b) unit cell lengths, c) intradimer dR,u_Ru dis- 
tance, and d) bond valence sums jl^ for Bi and Ru from 
Rietveld-refinement against NPD (Wombat@OPAL) data. 



by elongation of the stacking c axis (Fig. The 
transition appears to be first order in nature, with the 
slightly rounded disjunctions in refined structural param- 
eters being ascribed to the presence of both high- and 
low-temperature forms close to the transition (note that 
only a single phase could be refined at all temperatures, 
due to the resolution of the diffractometer used). In- 
terestingly, the intradimer distance rfRu-Ru — which lies 
along c — shows the opposite behaviour, shrinking from 
a maximum 2.76(3) A at 200 K to 2.66(3) A at 165 K. 

The apparent 4+ oxidation state for the single crys- 
tallographically unique Bi site in Ba3BiRu20g is very 
unusual (it would have an unstable 5d^^6s^ valence con- 
figuration), leading us to consider whether the structural 
transition is driven by disproportionation (2Bi^''" Bi'^^ 
-I- Bi^+) or charge transfer (Bi'^+ -I- 2Ru'^+ + 
2Ru^-^+). Either process should be obvious in bond va- 
lence sums (BVS) [Til calculated from experimentally 
refined structure parameters. BVS (Fig. [2Ji) show no 
anomaly for Ru on cooling, ruling out charge transfer. 
For Bi they shows a small anomaly; however, our NPD 
data show no evidence for symmetry lowering (additional 
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FIG. 3. (Color online) Magnetic susceptibility measurements. 
The solid lines are fits: eq. [5]- red line and eq. blue line. 
The bottom right inset shows the temperature dependence 
of the dx/dT derivative. The upper left inset displays the 
magnetic exchange interactions with nearest (Jo - within a 
dimer) and next-nearest (Ji, i = 1...5) neighbors. 



or split peaks below the transition) due to long-range 
Bi''+:Bi^"'' order. Moreover, Bi'*+ oxides generally dispro- 
portionate at much higher temperatures than the transi- 
tion of interest here {e.g., two distinct Bi sites can be dis- 
tinguished by BVS in BaBiOe even at 900 K The 6s 
electron on Bi''"'' probably localises at similarly high tem- 
peratures in Ba3BiRu209, but long-range Bi'^+:Bi^+ or- 
der is frustrated by the triangular disposition of Bi sites, 
with the only evidence being slightly anisotropic oxygen 
atomic displacement parameters (ADPs) (Fig. [T|). Most 
importantly, ADPs and diffraction peak widths show no 
discontinuities at the structural transition to indicate a 
significant change in local disorder or lattice strain. 

It does seem likely that "Bi'*+" plays an important role 
in the transition at T*, considering that no such transi- 
tion is observed for any i?*"*" {R — rare earth) cation. 
Unfortunately, the evidence at this stage is insufficient 
to identify the precise nature of that role. One possi- 
bility is that Bi^"*" may be acting as a transient charge 
reservoir, communicating changes in the electronic state 
between isolated RU2O9 dimers as they switch from the 
high- and low-temperature forms, thereby facilitating a 
first-order transition that never takes place where more 
stable i?'*"*' cations are involved. 

Molar magnetic susceptibility x is shown in Fig. [S] 
The structure consists of effectively isolated Ru20g units, 
so the expression for isolated S = 1 dimers: 



Xiso 



2jVAg'M|(l + 5exp(-2Jo/fcBr)) 
/cbT(3 + cxp(Jo/fcBT) + 5exp(-2Jo/fcBr)) 



(1) 



should be appropriate, where Na is Avogadro's num- 
ber, /iB is a Bohr magneton, fcB is the Boltzmann con- 
stant, g is the standard electron g ~ 2 factor and Jq is 
the intradimer magnetic coupling. However, a much bet- 
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ter quality of fit is obtained if interdimer interactions J' 
are included in the calculations in the mean-field approx- 
imation, so susceptibility Xdim can be written as: 



models for insulators, but the best fit was obtained with 



Xdi, 



Xu 



(1 + Axiso) 



(2) 



where A = J'/NAg^fJ-l- A fit above 180 K x(.T) = 
(1 — a;)xdim + Xo + xC/ (T — 9), which combines dimer, 
temperature independent and Curie- Weiss contributions 
(C = 1 emu K mole~^ for S* = 1 impurities) respectively, 
is shown. The best fit is to Jo/fce = 320(20) K, J'/fes = 
-160(10) K, a: = 0.02, 6 K and xo = 0.0011 emu 
mole~^. This is comparable to Ba3PrRu209f65, where 
Jo/k-Q =280 K. A comment should be made here con- 
cerning next-nearest neighbor interactions. While near- 
est neighbors are obviously those within RU2O9 units, 
higher-order neighbors are less obvious. Each Ru'*"'' is 
surrounded by 6 others in the pseudo-hexagonal ah plane 
and 3 from the neighboring plane, connected via Ru- 
0-Bi-O-Ru superexchange paths of approximately equal 
length. In the C2/c space group, these 9 sites com- 
prise 5 inequivalent positions, leading to 5 different ex- 
change constants Ji, i — 1...5 (inset in Fig. [3]). As 
J' = ''^i.Jii where is the multiplicity of next- 

nearest neighbors interacting via Ji, the average inter- 
dimer exchange parameter is Ji/k^ = —17.7 K. The in- 
tradimer coupling Jq is about double the absolute value 
of this interdimer J'. Susceptibility results thus show 
that above T* — 176 K, Ba3BiRu209 is an antiferro- 
magnetic (AFM) dimer system with weak ferromagnetic 
(FM) next-nearest neighbours interactions. 

A drastic drop in magnetic susceptibility below T* is 
highlighted by a maximum in the temperature derivative 
of x(T) (inset in Fig. [3]). As T* coincides with the dis- 
continuities in cell volume and cIru-Ru, Ba3BiRu209 ap- 
pears to undergo coincident structural/magnetic dimer- 
ization. We tentatively ascribe this to the opening 
of a gap in the spin-excitation spectrum between non- 
magnetic singlet ground state and excited triplet spin 
configurations within dimers. The magnetic susceptibil- 
ity of such a system should follow Q : 



Xsg = ar'^-'exp(— ^) 



(3) 



where a is a constant and A^^ is the value of the spin 
gap. A least-square fit performed in 2-150 K temperature 
range to the equation: x{T) = Xsg+Xo+c/(r— G) (shown 
in Fig. H blue line) yields a = 0.00021 emu mol-^K^^-^, 
A^/fcB = 247(2) K, Xo = 0.0006 emu mor^, c = 0.013 
emu mol~^K~'^, and Q = —5.5 K. 

Fig. H] shows the temperature dependence of electrical 
resistivity p. Ba3BiRu2 09 is nonmetallic - it exhibits 
increasing resistance with decrea sing temperature - like 
other Ru 6H-type perovskites [1, [l7f . We tested various 



the variable-range hopping expression 18j : 



los 



a + bT^ 



(4) 



where n is the dimensionality of hopping and a and b 
are constants. This transport mechanism is expected for 
an insulating sample with strong inhomogeneity, where 
charge carriers are localised into states with various en- 
ergies within a band gap. Transport between such states 
is realized with the help of a phonon, so that conduction 
occurs at a finite temperature. A least-square fit yields 
a = -5.2 rjcm, & = 99 K°-5f2cm and n = 0.97 [i.e. ~ 1), 
pointing to ID hopping. Below T*, p increases rapidly, 
seen as a drastic change in its temperature derivative (in- 
set in Fig. HI), and Eq. |4]is no longer valid. This is rem- 
iniscent of superconducting cuprates, where out-of-plane 
resistivity increases due to the opening of a spin-gap [l^ . 

Fig. [5] shows the temperature dependence of the 
heat capacity Cp/T ratio. The main feature is a sharp 
anomaly at T*, suggesting a large release of magnetic 
and/or structural entropy. The total heat capacity of a 
solid consists of lattice, magnetic, and electronic contri- 
butions, but without a suitable phononic reference mate- 
rial (isostructural but with no unpaired d electrons) the 
deconvolution of Cp is nontrivial. Initial attempts using 
a lattice term in the Debye approximation with a single 
characteristic frequency failed to produce a satisfactory 
fit. We have therefore used the approach of Junod et al. 
20l 21 1 , which has been successfull y appli ed to a number 
of oxide and intermetallic systems 22l424| , where the lat- 
tice contribution to the heat capacity is described by a 
sum of independent Debye and Einstein terms assigned to 
different groups of atoms. After testing all combinations 
of Debye and/or Einstein modes assigned to different sub- 
lattices, the best coefficient was obtained by fitting Cp 
below 150 K to a combination of 3 Debye modes assigned 
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FIG. 4. (Color online) Electrical resistivity on a semilog scale. 
The solid line is a fit to Eq. |4] The inset shows the tempera- 
ture derivative of p and its drastic change at T* = 176 K. 
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FIG. 5. (Color online) Temperature dependencies of the Cp/T 
ratio with its deconvolution into phononic and magnetic con- 
tributions (for details see the text). The inset shows apparent 
Einstein modes with e|" = 111 K. 



to the Ba, Bi and O sublattices, one Einstein mode from 
the Ru sublattice, and an additional term describing local 
magnetic excitations across the gap Ac 25|: 



Csg = nR 



(Ac/fcBr)^exp(-Ac/fcBr) 
1 -|-nexp(-Ac/fcBT) 



(5) 



where n = 3 for singlet-triplet excitations and R is 
the gas constant. We neglected the electronic contribu- 
tion because the sample was found to be semiconduct- 
ing. The best fit yielded eg' = 101(1), eg^ = 296(4), 
eg = 630(4), eg" = 111(1), and Ac/fce = 193(6) K. 
Note that Ac, estimated using this methodology, is of 
the same order as A^. These contributions are shown 
in Fig. [5l The inset shows the Cp/T^ ratio, where 
the maximum at 16 K is an apparent hallmark of Ein- 
stein modes [2^ 21 1. Neither a Debye model nor the 
spin-gap formula were able to reproduce this feature. 
The mean spin-gap from susceptibility and heat capacity 
A = (Ax + Ac)/2fcB w 220 K is equal to 0.7Jo/fcB. This 
is in perfect agreement with the A/ Jo ratio calculated 
for BaaMnzOs (A/Ajb = 12.2 K and JqAb =17.4 kQ), 
which has a similar structure consisting of hexagonal lay- 
ers of S' = 1 Mn5+ dimers. 

The transition at T* must be related to these changes 
in inter- and intradimer exchange interactions. We there- 
fore performed ab initio calculations on supercells for the 
structures refined at 300 and 2 K . As differences between 
total energies of different spin orientations appear due to 
the spin degrees of freedom, one can map total energies 
onto the Heisenberg hamiltonian: 



H — JijSiSj 

<«i> 



(6) 



with one intradimer exchange parameter Jq and two 
interdimer ones. Jab = (4Ji + 2J2)/6 (in-plane) and 



Jc = {J3 + J4 + Jd)/^ (along c) -please refer to the inset 
in fig. [3] At 300 K, intradimer coupling was estimated as 
J^^/k^ = 259 K in the presence of interdimer couplings 
jHT/^jg ^ _2.4 K and J^^ /k-Q = 5.4 K, pointing to the 
domination of exchange within RU2O9 units and its AFM 
character. At 2 K, the parameters are J}f^ /k-Q = 216 K, 
JLT/Zjg = 5.7 K and JlJ/k^ = 9.2 K. These values of 
Jo are somewhat smaller than obtained from susceptibil- 
ity (Jo/fcB = 320 K), but significant contributions from 
next-nearest neighbors {Jab and Jc) justify use of the in- 
teracting dimer model. It is worth noting that both AFM 
exchange integrals Jab and Jc at 2 K suggest magnetic 
frustration, which was proposed as the origin of the spin- 
gap opening in Mo3Sb7 |26l. |27|. 

In conclusion, we observe a spin-gap opening in the 
5 = 1 system BaaBiRuaOg below T* =176 K, with a 
spin- gap value of ~ 220 K. A magnetoelastic effect is 
observed as a decrease in the Ru— Ru distance within 
RU2O9 dimers below T*, as well as in the magnetic, ther- 
modynamic and electronic properties. An increase of unit 
cell volume is tentatively associated with a relaxation of 
the structure at T*, although we cannot exclude other 
possibilities such as orbital ordering. Ab initio calcula- 
tions confirm strong AFM coupling within RU2O9 dimers. 
Calculations to investigate the strong correlation effects 
typical of low-dimensional systems, using +Hubbard U 
methods, are now required, as well as direct confirmation 
of the spin-gapped state by inelastic neutron scattering. 

This work was supported by the ARC (DP0984585, 
DP0877695), AINSE, and the AMRFP. Dr Kevin Knight 
of ISIS assisted in HRPD data collection. 
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